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ABSTRACT 



We have calculated the one-loop supersymmetric corrections to the 2^2 parton scattering 
subprocesses qq qq, qq — > gg, gq — > gq and gg gg, including form-factor corrections and 
box diagrams with internal squarks and gluinos of arbitrary mass. In general, these exhibit 
cusps at the corresponding direct-channel sparticle thresholds. We use these calculations to 
make numerical estimates of the possible threshold effects at the Fermilab Tevatron collider 
and at the LHC, which depend on the rapidity range selected, but can be as large as a few 
percent. These effects are diluted in the integrated large- ii^-r cross section, where they are 
negative. 
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1 Introduction 



Large samples of inclusive liadronic jets are available now from the Fermilab Tevatron col- 
lider, and will become available from the LHC in the future. These offer the prospects of 
precision tests of QCD, which is of interest and importance in its own right. Moreover, as 
the dynamics of QCD is better understood, it becomes a sharper tool for probing possible 
physics beyond the Standard Model. Complete calculations of jet physics to next-to-leading 
order in perturbation theory exist in pure QCD . However, this may not be sufficient for 
precision physics at an exploratory machine that opens up a new energy range where there 
may be thresholds for new physics. For instance, the Fermilab Tevatron collider has already 
crossed the threshold for it production, and either it or the LHC should cross the threshold 
for squark and gluino production. In the neighbourhood of, and beyond, the threshold for 
such new heavy strongly-interacting particles, their virtual effects should be included in a 
complete treatment of one-loop perturbative QCD effects. 

Our attention was drawn to this problem by recent measurements of large- jet cross 
sections at the Tevatron collider 0, ^ , some of which exhibited a prima facie discrepancy with 
predictions based on the parton distribution functions available previously. Clearly there are 
uncertainties in these distributions P], which are not (yet) calculable from first principles, 
and it has been argued in particular that reasonable uncertainties in the gluon distribution 
could accommodate simultaneously both the CDF and DO jet measurements. However, 
it has also been argued that strong-interaction threshold effects could be significant ^ |^, 
with emphasis placed on the possible importance of the sparticle threshold |§] . One might 
consider whether such a threshold effect could provide an alternative signature for sparticles, 
if their decays differ from those anticipated in conventional direct searches. Even if one does 
not expect the sparticle threshold to be very important, it is clearly desirable to have an 
exact one-loop treatment of it, just as one- loop sparticle corrections are known and used in 
the analysis of precision electroweak physics at and around the Z° peak 

We have undertaken detailed calculations of such corrections in the Minimal Supersym- 
metric extension of the Standard Model (MSSM), and published first results [0. We found 
that there could be threshold structures of the order of a few % in the cross sections for 
individual parton-parton scattering subprocesses when one-loop corrections were included. 
However, these were of the wrong shape and of insufficient magnitude to make a significant 
contribution to the resolution of the CDF large- conundrum 0. The purpose of this paper 
is to complete the analysis of ||T^, presenting more details of the calculations, presenting the 
one-loop corrections to all partonic subprocesses, including box diagrams and the gg —>■ gg 
subprocess, which were not included in and combining our calculations in a complete 
one-loop numerical analysis Q of jet cross sections at Fermilab and the LHC, including the 
appropriate convolutions over all the contributing parton-parton scattering subprocesses. 

The layout of this paper is as follows. In section 2 we present details of our results |0 
for the calculations of the one-loop corrections to the partonic subprocesses arising from the 
form-factor corrections to the triple-gluon vertex and the quark-quark-gluon vertex, which 
include self-energy corrections due to squarks and gluinos. We comment on various consis- 

^ Numerical results of such a study have also been published in [pl| . 

^We assume for simplicity that all the relevant left- and right-handed squarks have equal masses. This 
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tency checks on our calculations, demonstrating in particular cancellations between different 
diagrams associated with Ward identities, and remaining logarithmic singularities that reffect 
the expected running of the strong coupling as above the supersymmetric threshold. Section 
3 contains a similar analysis of one-loop box diagrams, including a number of further consis- 
tency checks. Numerical results are presented in section 4, including the convolutions with 
parton distribution functions which are appropriate for the analysis of jet cross sections at 
the Fermilab Tevatron collider when ~ m = = mg = 200 GeV (where rUs^g denote the 
squark and gluino masses), and at the LHC when Ex ~ m = 1 TeV. Section 5 summarizes 



our conclusions from our calculations. In the appendix we list the Veltman-Passarino ^2 
functions, in terms of which our results are displayed. 



2 Form-Factor Corrections 



In this section present the form-factor corrections for the partonic subprocesses 

q{pi) + q{p2) qips) + q{pi) (2.1) 
q{pi) + q{p2) g{pz) + 9{Pa) (2.2) 
9{Pi) + 9{P2) g{pz) + 9{Pi)- (2.3) 
calculated using the MS prescription. The form factors for all other partonic subprocesses 
can be obtained from these by exploiting crossing symmetry. 

The contribution to the amplitude for the subprocess ( p.l|) from all diagrams except box 
graphs may be written in the general form 



viP4)F,,gis)rrMP3)viP2)F'^'^'is)rr''uiP 



1^ 



(2.4) 



where is the generator of colour SU{3) in the defining representation, and Fggg{s) is a 
form factor for the quark-quark-virtual-gluon vertex. The tree-level amplitude is obtained 
by replacing this form factor by unity. The one-loop corrections to this form factor from 
internal squarks and gluinos come from the diagrams shown in Fig. [I|(a,b,g,h), and give a 
contribution to Fggg of the form 



AF, 



Qqg 
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C24({1}) + Ca 



+ 2C24({2})-m2Co({2}) 



+s {C,,{{2}) + Cum)) - - 



1 /2m' 



— + lUBo{{l}) + l) + -A{mg) 
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^ (5o({2}) + l) + -A(m, 
SI s 



(2.5) 



where the Veltman-Passarino (VP) functions [T^ Cj, Bq are given in the appendix. They 
are calculated in n = 4 — 2e dimensions, and the arguments {1}, {2} are defined by 



Ci({l}) = Ci{s,ms,mg,ms), Ci{{2}) = Ci{s,mg,ms,mg), 
Bo{{l}) = Bo{s,mg,mg), Bo{{2}) = Bo{s,ms,ms), 



(2.6) 



means that our results should not be assumed to apply to processes involving both sbottom and stop squarks. 
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(i) 



Figure 1: One- loop Feynman diagrams involving virtual sparticles in the MSSM for (a,b) 
the viitual-gqq vertex, (c,d) the yiitual-ggg vertex, (e,f) the viTtual-qqg vertex, (g,h) the 
gluon self energy, and (i) the quark self energy. Here and in subsequent figures of Feynman 
diagrams, the broken lines represent squarks and the double solid lines represent gluinos. 
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where and nig are the squark and gluino masses, respectively. For the unsubtracted 
correction to the form factor, we must add to (|2.5| ) the pole term 



- C^-, (2.7) 
e 



which arises from the wave-function renormalization of the external quarks, and cancels the 
pole part of ( p.5|) proportional to Cp, as required by the abelian Ward identity. Substituting 
the ultraviolet-divergent functions by their pole parts, given in ( |A.l^j| - |X?T3| ), we see that the 



poles from the one-particle irreducible vertex corrections shown in Figs. |l|(a,b)) which are 
proportional to Ca also cancel between the two diagrams. The the remaining divergences 
come from the functions i?o, which contribute to the extra renormalization of the strong 
coupling due to supersymmetric particles, namely 

{Ca + Tr), (2.8) 
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and arise from the gluon self-energy contributions shown in Figs. |I](g,h)). 

Once this renormalization has been effected, the finite form factor is given by 
with the divergent functions Bq and C24 understood to have been subtracted using the MS 
prescription. The resulting contribution to the form factor vanishes in the limit s ^ when 
we set TUs = mg = m. In the limit s <^ m? the form-factor correction is given by the 
following simple expression: 



,,..9 119 /i„ 1 o «n on ™2 ^ ' 



/<<m2 ~^ ''''^ 477 I 12 60 30 J m2 ' 

which is 0(O.l(s/m^))%. This and other form-factor corrections are actually numerically 
considerably larger when s ~ [|l^ , as we shall see later. 

The amplitude for the process ( p.2|) , again omitting for the moment box diagrams, may 
be written as 

- lAT^a^ey, |^t;(p2)r"^r'^^F;,,(pi - VzWi - h )i^;,,(Pi - ^3)^(^1) 

-lv{p2) [r^^r^^] Fg,g{s)YVr{P3,P4Hp,)^, (2.10) 

where €3, €4 are the polarization vectors of the external gluons with momenta p^, p^, whose 
colours are 03, 04, respectively. 

The quantity F^g^{q) is the form factor for the gluon-quark- virtual-quark vertex, which at 
the tree level is simply 7^^. The diagrams contributing to the corrections to this form factor 
which involve internal squarks and gluinos are shown in Fig. |l|(e,f,i), and give a contribution 

^KM = ^ii^cF-^CA)[C2,mh' + iC22m)-c2smw] 



47r 
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+Ca 



2C24({5}) - - - mJCo({5}) - q\CM^}) + CMh})) ) + 2CMm<l' 
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5o(g ,ms,m3) - — (A(mg) - y4(m^)) 

where the arguments {4}, {5} are given by 

Cj({4}) = Ci(g^ m,, m^, m^), Ci({5}) = m^, m^, mj. 

To obtain the unsubtracted form factor, we must add a pole term 

as 



r 



Aire 



Cf Ca Tr 

T ^ T 



f2.111 



(2.12) 



(2.13) 



which arises from the wave-function renormahzations of the external quark and gluon legs. 
Once this is included, we notice that again the abelian pole term proportional to C^? cancels. 
The remaining pole term is given by ( |2.8| ), and is absorbed by the renormalization of the 
strong coupling. After this renormalization, and setting = = m, we find that for 

1 , 



q^ 
Am? 



+ 



(2.14) 



which is (9(l(g^/m^ 



The three-gluon vertex function V'^^''{p'i,pA) may be written as 



(2.15) 



where s = 2p3 -p^. At the tree level, the form factors Fi, F2 take the value unity, whereas F3 
is zero. The irreducible one-loop contributions from squark and gluino loops to these form 
factors are are shown in Fig. |l|(c,d), and are given by 



AFi(s) 



Att 



-3-K(^o({6})-Cn({6}) + Ci,({6})) 



+ S (Ci2({6}) + C22({6}) - C33({6}) + C34({6})) 
-2C35({6}) + 2C36({6}) + 2C24({6}) 

1 (2ml ,\ 2 . A 

~3 



-^ + 1 (i?o({l}) + l) + ^A(m,) + - 



+TJ 



R 



1 / 4m^\ 4 
8C36({7}) - 8C35({7}) - - 1 - — i (5o({2}) + 1) - -Aim, 



(2.16) 



AF2(s) 



47r 



-3-K(Co({6})-Ci2({6})) 



-S (C22({6}) + C34({6})) - 2C36({6}) + 2C24({6}) 

1 (2ml .\ 2 , A' 

'3 



-l + lUB,{{l}) + l) + fA{m,) + - 



[B,{{2}) + 1) - ^^K) + g 



(2.17) 
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AF3(.) = ^|C^[C22({6})-C23({6})-C33({6})+C34({6})] 



TT 



+2Tr [C23({7}) - C22({7}) + C33({7}) - C34({7})] |, (2.18) 

where the arguments {6}, {7} are given by 

= Ci{s, Trig, nig, nig), Ci{{7}) = Ci{s, rris, rris, m^). (2.19) 
To obtain the unsubtracted form factor, we must add a pole term 



/2Ca , 2T, 



+ =^ ) , (2.20) 
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which arises from the wave-function renormahzations of the external gluon legs as shown in 
Fig. |l|(g,h), leaving a pole term equal to ( p.8|) . After renormalization and setting = rrig = 
m, we find that for s <^ 

lim AFUs) ^ — (2.21) 
hm^ AF2(.) - ^ (Ca + ^) 4 (2-22) 

(T 1 

hm AF3(.) ^ (C^ - 2Tn) — (2.23) 
which are comparable in magnitude to the low-energy expansions of the previous form factors. 



Apart from some minor differences, the low-energy expansions of our results (eqs. (|2.9| , p.l4| . 



|2.21| , ^.22| , |2.23| )) confirm the general magnitudes of the below-threshold corrections found 



m ^. 

Finally, we can use V3 to write down the amplitude for process ( p^.3D excluding the 
box diagrams which contribute to the renormalization of the four-gluon coupling, which 
we postpone until the next section: 

,l!L^e^^6geI/"^'^^V'^^"^'V^3''''^(pi,P2)V3^'"^(p3,P4) + permutations {2,3,4}, (2.24) 

s 

where ei, €2 are the polarization vectors of the incoming gluons with momenta pi, p2, whose 
colours are ai, 02, respectively, and €3, are the polarization vectors of the outgoing gluons 
with momenta p3, ^4, whose colours are 03, 04, respectively. The notation "permutations 
{2,3,4}" means permutations over the momenta, polarization vectors and colours of gluons 
2,3,4. 



3 Box Diagrams 



For the box diagrams we need the four-point VP functions Di [|T2] listed in the appendix in 
O to ( [AH . 
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Figure 2: Box diagrams for the process qq ^ qq. 
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Figure 3: Box diagrams for the process qq 
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(c) 



mum^P't 
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The box diagrams for the process ( |2.1| ), shown in Fig. give the following contribution 
to the amplitude: 

2iCAal {D2{t, u, rris, rUg, rris, rUg, fx, v) [A'^y + A''^^ 

+mlDQ{t,u,ms,mg,ms,mg)[Ass + App]^ + [t ^ u) (3.1) 

where 

= v{p,)i^T%{p2) u{P3)Yr^uip,) 
A7a = viP4)Yl'r%iP2) uiPsh'l'r'^uip,) 
Ass = v{p4)t''v{p2) u{p3)t''u{pi) 

App = v{p4h^T%{p2) u{p3)-f^T''u{pi) 

The second graph in Fig. |^ arises from the Majorana nature of the internal gluino lines, 
which permits them to propagate both along and against the flow of the quark fermion 
number. 

The box diagrams for the process ( |2.2|) are shown in Fig. ^ Fig. ^(a) gives a contribution 

- 2^0^6364 { [ D-sit, u, nis, rrig, nig, nig, /i, p, r) 

+ {Pi - P3yD2it, u, nis, nig, nig, nig, /i, r)] B^l^^'''' 
+ nil[Di it, u, m„ nig, nig, nig, /i) {B^ + + K7) 

+ (Pi - PsTDoit, u, nis, rrig, nig, mg)B'X''] } +(3^4), 
where the tensors B are defined as 



r)fj.l--fj.n 



T ,r 



r ,r 



u{pi 



■V{j)2) (7'"^ ■■■7- 

Fig. |^(b) gives a contribution 

- Sialeltl { D^{t, u, nig, m^, nis, rUs, fi, v, a) + plD2{t, u, nig, nis, rris, nis, cr) 
+ {Pi -P^Y [D2(t,u,mg,ms,ms,nis, fi, v) + p\Di{t,u,mg,nis,ms,nis, ^i)]} 
v{p2) (7'^rVVV^) u{pi) +(3^4) 



(3.2) 
(3.3) 



(3.4) 
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Figure 4: Gluino box diagram for the process g g ^ g g- 
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Figure 5: (a) Squark box diagram and (b,c,d) diagrams involving quartic couplings that 
contribute to the process gg ^ g g- 

Fig. |^(c) gives a contribution 



- Aialelel { [ D'i{t, s, rris, nig, nig, m^, /i, p, cr) 
+ (Pi - P3TD2{t, s, rus, rUg, rUg, rUs, p, a)] B'^^"'' 

+ mlDi{t,s,ms,mg,mg,ms,a)B'^^ } + (3 ^ 4), 

where the tensors B' are defined by 

and the notation (3 ^ 4) means (es <-> €4, P3 ^ P4, t ^ u, a ^ b). 



(3.5) 



(3.6) 



The contributions to the process ( [^.3D from box diagrams are shown in Fig. ^ (for an 
internal gluino loop) and Fig. ^ (for an internal squark loop). Also included in Fig. ^ are 
diagrams that involve the four-point coupling between gluons and squarks. Although these 
are not box diagrams, we include them here, because they do not contribute to the form 
factor for the triple-gluon vertex. The contribution to the amplitude from Fig. ^ is 

Gl { Di{s, t, rrig, rrig, nig, nig, /i, u, p, cr)tT (^ 7^3 7^") 
+03(3, t, nig, nig, nig, nig, p, V, p) [tr (^ Yh^^k Y) 

+tr(A7'^A773 7"^^i) +tr(A774 773(^i +^2)^27^ ] 
+1^2(5, t, nig, nig, nig, nig, p,u) [tr (^1 7^ A ^4 ^ 7^"^ l^i ) 
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+Di{s,t,mg,mg,mg,mg, iJ,)tT {^i'y''^i:^i^s (^i +^2)^^) 
+mlD2{s, t, nig, nig, nig, nig, /i, z/) [16e'(e'^e2 ■ £4 + We^el^ei ■ eg] 
+nilDi{s,t,nig,nig,nig,mg,fi) [2e^tr(^^ (^1 + ^2)^) 

+2e^tr (/^i ^ ^4 ^ ) + 2e^tr (^1 ^ ^ ) + tr i>A h h ) 

+tr ^ ^ 1^1 ) + tr (77i ^ ^ + ^2 ) ) ] 

+mjL'o(s,)f:,mg,mg,mg,mg) [tr (^^^4^ (^1 + ^2)^) 
+tr ^ ^ (^1 + 1*2 ) ^ 1*1 ) + tr 64/4 ^ ^ ^1 ) ] 
+m^L'o(-s, t, nig, nig, nig, mg)tr ^ ^ ^ ) } 

(+ permutations {2, 3, 4}), (3.7) 



In the interests of compactness, we have left this expression in terms of traces over 7 matrices. 



where the colour factor Gi is given by 



j^aiab ja2hc j^a-^cd ja^da 



(3.8) 



The contribution from Fig. ^(a) is 



32ia^G'2 (^1(^24(^1 { Da{s, t, m^, m^, m^, m^, /i, z/, p, a) 



+p'(D3{s, t, nis, nis, nis, m,, /i, p, a) 

+ (Pi + P2Y Ds^s, t, nis, nis, nis, m^, p, z/, cr) 

PiiPi + P2YD2{s, t, nis, rris, nis, rUs, p, a) } 



+ permutations {2,3,4}), 



(3.9) 



where the colour factor G2 is given by 



G2 = tT{T''W^T''^T''^) 



(3.10) 



The contribution from Fig. ^(b) is 



Aia^Gs ei ■ £4 €2 ■ €3 Bq{u, nis, nis) 
(+ permutations {2,3,4}), 



(3.11) 
(3.12) 



where the colour factor G3 is given by 




(3.13) 



The contribution from Fig, ^(c) is 

-IGia^Gs 62 ■ ea [ C24.{u, nis, rris, nis) ei ■ €4 - C2z{u, nis, rris, nis) Pa • ei pi • €4] 
(+ permutations {2,3,4}), 

Finally, from Fig. ^(d) we get 

-IGia^Gaei ■ €4 [C2A{u,ms,nis,nis) €2 ■ €3 - G2z{u,nis,nis,nis)pi. ■ €2^2 ■ £3] 
(+ permutations {2,3,4}), 



(3.14) 



(3.15) 
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These contributions contain ultraviolet divergences, which are associated with the renormal- 
ization Z4 of the four-point gluon vertex. Application of Ward identities tells us that the 
coefficient of the pole parts of the integrals should be given by the contribution A(3q to the (3 
function, multiplied by the tree-level amplitude r4 for the four-gluon coupling. In the above 
expressions, poles arise in the unsubtracted forms of the functions Bq, C24-D4. Exploiting 
the pole parts displayed in the appendix in ( |A. 14|JX?T3| , |A. ISj ) , and permuting the gluons 
2,3,4, we arrive at a pole term 

n 2 1 

- {CA + Tn)T, = -A(3oT^, (3.16) 



47re3 ' ' e 

as required. 

In order to obtain the relevant contributions to the partonic cross sections, these contri- 
butions to amplitudes must be multiplied by the Hermitian conjugates of the corresponding 
tree-level amplitudes and summed (averaged) over final (initial) quark or gluon polariza- 
tions and colours. These tedious but straightforward manipulations are most conveniently 
performed using a fast algebraic manipulation package: we have used FORM. 

Finally, we note that the box-diagram contributions for all other partonic subprocesses 
can be obtained from the above expressions by crossing symmetry. 



4 Numerical Results 

We first present the results for the one-loop corrections to all the different parton-parton 
scattering cross sections shown in Fig. These are for the simplified case = nig = m, 
although our analytical results apply for arbitrary ratios rris/mg. These extend and complete 
the analogous plots shown in [|l^], which did not include box diagrams, nor results for the 



subprocess gg — >• gg. We already commented in |T0[ that we did not expect the box- 
diagram contributions to be numerically large, and this has been confirmed by our complete 
calculation. 

We first note that, because these corrections to the cross sections are due to interferences 
between tree-level and one-loop amplitudes, they are not necessarily positive. In particular, 
processes which are dominated by the exchange of a parton in the t channel give a negative 
correction, and we see that the corrections in Fig. P(c,d,e) are negative for s ~ 4m^. We also 
note that each of the corrections exhibits a cusp at ^/s = 2m, creating a local maximum in 
the magnitude of the one-loop correction. This is largest for the subprocesses qjqj — > qkQk 
and qq gg shown in Figs. ^(a,b), which are unfortunately not dominant at the Fermilab 
Tevatron collider and the LHC Next in magnitude is the qjg qjg shown in Fig. which 
is of greater experimental significance. The newly-calculated one-loop correction to the 
subprocess gg — > gg shown in Fig. ^(e) is the smallest numerically. 

We also note that each of the one-loop corrections grows logarithmically at large s ^ m^. 
As was mentioned previously, some of this logarithmic behaviour is due to ultraviolet diver- 
gences associated with the different running of a<j above the sparticle threshold. However, 
this is not the only source of logarithmic behaviour. Since the two-jet cross section does not 
include final states containing sparticles, there are also infrared divergences when -C s. 
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Figure 6: One-loop virtual-sparticle corrections in the threshold region of the subprocess 
centre-of-mass energy squared s to the processes (a) qj qj — > g^, (b) qq ^ g g, (c) qj q^ — > 
qjqk, (d) qj g qjg, and (e) gg — >■ gg. In each case we have plotted the ratio of the 
one-loop cross section to the tree-level value, and in (b,c,d,e) we have plotted this ratio for 
three different choices of t/s, equivalent to different parton-parton centre-of-mass scattering 
angles. All corrections are evaluated using as = 0.11, and = rUg = m. 
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Moreover, numerical studies show that the logarithmic asymptotic behaviour does not set in 
until subenergies beyond the reach of the Fermilab Tevatron collider - for m = 200 GeV - or 
the LHC - for m = 1 TeV. For these two reasons, it is not adequate to model the sparticle 
threshold simply by switching to the above-threshold form of for s > Am? . 

We now use the above results to calculate one-loop corrections to jet cross sections at 
large Ex and two-jet invariant mass M, by convoluting the above subprocess results with 
parton distributions. Since the higher-order sparticle corrections are only significant over a 
small range of parton subenergy in the threshold region, it is convenient to start with the 
triply-differential cross section as a function of M and the rapidities ?/i,?/2 of the final-state 
jet pair. In terms of the partonic squared matrix elements, this is given by 

'^''^ ^ E/^(^i)/.(^2) A<g(M^^)p + |>^g(M^^)r) (4.1) 



ijkl 



where fi{x) is the parton distribution function for parton i, Ai^j is the matrix element for 
the scattering of parton i and parton j to partons k, I, and the jet rapidities yi^2 are given 
by 

x^2) = ^e^^y^^y-\ (4.2) 

and 

t=-M'-u = -— -. (4.3) 

1 +exp(?/i -y2) 



To exemplify the results obtained with (|4.1| ), we consider a number of discrete choices of the 



jet rapidities ?/i,2, shown in Fig. ^for the Fermilab Tevatron, and in Fig. ^for the LHC. We 



have used the parton-distribution functions from ||13[, and do not expect that using other 
parametrizations would significantly affect our results. 

The first striking feature is that the corrections are negative. This is because, over most 
of the ranges of parton momentum fractions Xi, X2 studied, the subprocesses involving t- 
channel exchange dominate over the annihilation processes. Indeed, recalling that we must 
symmetrize over the partons in the final-state jets, we see from the tree-level cross-sections 
given in that the subprocess cross section for the scattering of two quarks with different 
flavours at any value of s is always an order of magnitude larger than the corresponding 
cross section for quark-antiquark annihilation, quite apart from the fact that the parton 
distribution functions provide more flux for this process. The second important feature of the 
plots shown in Figs. 05§|is that, although characteristic cusps always appear at threshold, the 
reduction of the differential cross-section due to the sparticle-loop corrections is significant 
only for narrow ranges of rapidity pairs, {?/i,?/2}- For this reason, these corrections are 
substantially washed out if one integrates over one of the rapidities, in order to obtain the 
double-differential cross-section, in which the rapidity of only one of the final-state jets is 
measured. Nevertheless, for sufficiently small values of jet rapidity, there still appears a cusp 
at threshold, as shown in Fig. ^ for the case of the Fermilab Tevatron collider. However, 
these cusps are always somewhat diminished and broader than the previous cusps for the 
triple-differential cross sections. This "washing out" effect is greater for the LHC shown in 
Fig. ITOl, where the dip in the threshold region is quite possibly too small to be observed. 
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Figure 7: Triple-differential cross section for pp scattering at ^ =1.8 TeV, as a function 
of the two-jet invariant mass M, for various choices of jet rapidities, calculated assuming 
"i^s = fT^g = 200 GeV. Here and in subsequent plots, R denotes the ratio of the one-sparticle- 
loop-corrected cross section to the tree-level cross section. 
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Figure 8: Triple-differential cross section for pp scattering at y/s — 14 TeV, as a function 
of the two-jet invariant mass M, for various choices of jet rapidities, calculated assuming 
nis = rUg = 1 TeV. Again, R denotes the ratio of the one-sparticle-loop-corrected cross 
section to the tree-level cross section. 
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Figure 9: Double-differential cross-section for pp scattering at =1.8 TeV, as a function 
of the two-jet invariant mass M, for various different jet rapidities, calculated assuming 



= = 200 GeV. 
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Figure 10: Double-differential cross-section for pp scattering at ^ — 14 TeV, as a function 
of the two-jet invariant mass M, for various different jet rapidities, calculated assuming 

T^s = TTT'g = 1 TeV. 
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Figure 11: Single-differential cross section for pp scattering at ^/s = 1.8 TeV, as a function 



of transverse energy Et-, calculated assuming 



200 GeV. 



This dilution is even more evident in the case of the single-differential cross section as a 



function of the transverse energy Et^ as shown in Fig. ITT] for the case of the Tevatron. There 
is a broad dip in the cross section for Et — m, but this is probably also too shallow to be 



observable. This dip is again shallower for the LHC large- E't cross section, shown in Fig ^ 
In general, the sparticle-loop signal is clearer for sparticles of mass ~ 200 GeV at the 
Tevatron than for sparticles of mass ~ 1 TeV at the LHC. The reason for this is that smaller 
values of m^/s are sensitive, on the average, to smaller values of x, thereby sampling more 
of the gluon content of the incident hadrons. As we saw earlier, the sparticle-loop effects 
are smallest for purely gluonic scattering. It is, however, noteworthy that the net effect of 
virtual-sparticle loops is to decrease the predicted cross section, and can therefore not be 
used even as a partial explanation of any unexpected rise in the large- differential cross 
section 0. 



5 Conclusions 

We have presented in this paper complete one-sparticle-loop corrections to the large- i^T and 
large-M cross sections at high-energy hadron-hadron colliders, and used them to present 
some numerical results for the Fermilab Tevatron collider and the LHC. We find that the 
sparticle-loop effects are too small, and of the wrong sign, to make any contribution to 
explaining the possible laige-Ex cross-section discrepancy reported recently 0. However, 
the fact that our calculated corrections exhibit cusps at the sparticle threshold may encourage 
the hope that these effects could be visible in the large-statistics data to be obtained in the 
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Figure 12: Single-differential cross section for pp scattering at ^/s — 14 TeV, as a function 
of transverse energy Et, calculated assuming = nig — 1 TeV. 

future at the Tevatron collider and the LHC. 

As we have shown, these cusps are most noticeable in the triple-differential jet cross 
section d^o'/dMdyidy2, and get progressively more washed out as one integrates over one or 
both jet rapidities, or if one plots the integrated large- £'r cross section. It is for the experi- 
mental collaborations to judge whether they will be able to obtain the necessary statistics, 
and whether the systematic errors can be controlled to the desired low level. In this paper 
we have not included any allowance for experimental effects such as the initial transverse 
momenta of the colliding partons, extra gluon radiation, or the experimental resolution in 
the large- ii^T jet energies. 

Evidently, we do not know where the squark and gluino thresholds may be, nor whether 
they are coincident. The numerical results presented in this paper have been for the op- 
timistic case nis = mg, and our threshold effects will be surely be spread out and further 
diluted to some extent if and rrig are substantially different. However, we would like to 
point out that looking for such a cusps in differential cross sections is in principle a model- 
independent, though indirect, way of looking for strongly-interacting sparticles. The only 
vertices that enter our calculations are those proportional to ctg, and the squark and gluino 
decay vertices do not enter. Thus, looking for the subtle effects we have presented here is 
a strategy complementary to the direct searches for sparticles decaying according to some 
particular model scenario, which depends whether R parity conserved or not, and on the 
spectrum of lighter sparticles. Also, the large cross section for large- jets provides a win- 
dow on large-mass physics that may reach out to larger sparticle masses than direct searches 
in specific decay modes with uncertain branching ratios. 

Anybody interested in obtaining the code used to derive these results should contact 
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D.A.R., who will make it freely available. 
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Appendix 

In this appendix we list the VP |jl2| functions used in the text. 
Tadpole function 

Two-point Functions: 



Vertex functions: 

d^k 



Co (A.4) 



(27r)4 (p -mf)((fc + pi)2 -m2)2((fc + pi +P2)^ - mi) IQn 

f ^ ^ J- ^ 

J (27r)4 (P - ml){{k + pi)2 - m2f{{k + pi+ P2Y - mfj ~ IQti^ ^ 

(A.5) 



d'^k ki^k" 



{271^ (A;2 - ml){{k + pif - m2f{{k + pi+ P2Y - mi) 

T^2 (^21 + C^22P^P^ + C2MV2+f2P\) + ^24^7^''') , (A.6) 



i (27r)" - m2)((A: + p^f - m2Y{{k + pi + ^2)' - ml) 

( C31 + C32 P^P^P^ + C33 (p + p^'p^p? + pM) 

+C34(PlP2P2 + + P2P2P1) 

+ C^^{g^^p{ + ^^pPi + ^^pp^*) + C3e{9,^uP2 + 9fipP2 + ^!^/#2) ) , (A.7) 

where the functions Co, Cij have in general the arguments (Pi,P2, (pi + P2Y ,m\,m\,m\). 
However, in our case we always have Pi = P2 = 0, so we suppress the first two arguments. 
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Box functions: 



dt^k 



Do{2pi ■P2,2p2 ■P4,mi,m2,m3,m4), 



ki" 



(A. 



(27r)'^ (/c^ — ml){{k + piY — mDHk + pi + p-if — vn^[{k — ^4)^ — m\ 



167r2 



d^k 



Di{2pi ■ P2, 2p2 ■ Pi, rni, ms, ms, rrii, /i), 



(A.9) 



(27r)'^ _ rnl){{k + pi)^ — m2)((/c + Pi + ^2)^ — "^Dll^ ^ ^4)^ ^ "^1 



167r2 



^2(2pi ■ P2, 2p2 ■ Pi, rni, m2, m^, rrii, /i, u) 



k^'k'^kP 



(A.IO) 



(27r)4 (P - mf)((A; + - m^)((A; + pi +^2)^ - - ^4)^ - "^1) 

i 



167r2 



-D3(2pi ■ p2, 2p2 ■ Pi, mi, m2, ma, m4, /i, i^, p), (A.ll) 



r d^'k 

J f27r)" - 



k/^k^kPk" 



(27r)"- (/c^ — rnl){{k + piY — mDdk + Pi + P2f — vn^[{k — ^4)^ — ml) 



;D^{2px ■ p2, 2p2 ■ Pi, mi, m2, ms, m4, /i, z/, p, a) 



(A.12) 



For the sake of compactness, we do not write these tensor expressions out in terms of the 
vectors pi ■ ■ - pi, but refer the reader to [|12| for details. 



The exact forms of the functions Bi,Ci, Cij, Di are given in ||12[ . The functions A, Bq, Bi, C24, 
C35, C36, -D4 are ultraviolet divergent, and therefore should be calulated in n = 4 — 2e dimen- 
sions. The pole parts of these functions are given by the following expressions: 

i2 



m 



P.P{A(m)} = — 



1 



P.P{5o(x,mi,m2)} ^ 

P.P {C24(a;, mi, m2, mg)} = 

P.P {C35(x, mi, 1712, ms)} = - 

P.P {Cseix, mi, m2, mg)} = - 
1 



1 

1 

1 



P.P {^4(0:, y, mi, m2, m^, m,, p, iy,p,a)} = — (g'^-'g"'^ + g'^'^g'"^ + g^'^g-'^) 



(A. 13) 
(A. 14) 
(A.15) 
(A.16) 
(A.17) 
(A.18) 
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